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We introduce the heterogeneous-A;-core, which generalizes the fc-core, and contrast it with boot- 
strap percolation. Vertices have a threshold ki which may be different at each vertex. If a vertex 
has less than ki neighbors it is pruned from the network. The heterogeneous-fc-core is the sub-graph 
remaining after no further vertices can be pruned. If the thresholds ki are 1 with probability / or 
fc > 3 with probability (1 — /), the process forms one branch of an activation-pruning process which 
demonstrates hysteresis. The other branch is formed by ordinary bootstrap percolation. We show 
that there are two types of transitions in this heterogeneous-fe-core process: the giant heterogeneous- 
fc-core may appear with a continuous transition and there may be a second, discontinuous, hybrid 
transition. We compare critical phenomena, critical clusters and avalanches at the heterogeneous-fc- 
core and bootstrap percolation transitions. We also show that network structure has a crucial effect 
on these processes, with the giant heterogeneous-fc-core appearing immediately at a finite value for 
any / > when the degree distribution tends to a power law P{q) ~ q~~' with 7 < 3. 
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Bootstrap percolation and the fc-core are closely re- 
lated concepts, and in fact it is easy to confuse the two. 
Both belong to a new class of systems with hybrid phase 
transitions, yet it can be clearly shown that the two pro- 
cesses do not map onto each other. Here we elucidate 
the relationship and differences between these two con- 
cepts by introducing a generalization of the /c-core, the 
heterogeneous-fc-core. 

The fc-core is the maximal sub-graph whose vertices all 
have internal degree at least fc It has proved a use- 
ful tool giving insight into the deep structure of complex 
networks [2rM ; ^^'^ found applications in diverse 
areas, from rigidity 01 and jamming Q transitions to 
real neural networks [3, and evolution flU . The k- 
core has been extensively studied on tree-like networks, 
starting with Bethe lattices [H, [l^ and Random graphs 
before finall y b eing extended to arbitrary degree 
distributions [1, IT7l4l9| |. Hyperbolic lattices have also 
been considered [20|. Other studies, mostly numerical, 
have considered the sizes of culling avalanches [2l| - [23j . 
Results on non tree-like graphs have been largely numer- 
ical [13, HE] ) although some analytic results incorporat- 
ing clustering have recently been obtained [1^ [23| • At 
the same time, bootstrap percolation has emerged as a 
useful model for a variety of applications such as neu- 
ronal activity [28l - [30|. ja mming and rigidity transitions 
and glassy dynamics [3l|, [s^ , and magnetic systems [s^l ■ 
In bootstrap percolation, a set of seed vertices is initially 
activated, and other vertices become active if they have k 
active neighbors. This process has been investigated on 
two and three dimensional lattices (see [s^ - fSTj and ref- 
erences therein). Bootstrap percolation has been studied 
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on the random regular graph [38|, |39j, on infinite trees 
poj . and most recently general complex networks [4l| . 
Finite random graphs have also been studied [13]. An 
interesting alternative formulation is the Watts model of 
opinions, in which the threshold is defined as a certain 
fraction of the neighbors rather than an absolute number 
[4^ . These processes may also be generalized so that the 
thresholds may be different at each vertex [U . 

Here we introduce a generalization of the fc-core, the 
heterogeneous-fc-core. In the heterogeneous-fc-core, each 
vertex i in a network has a hidden variable, its threshold 
value ki. The heterogeneous-fc-core is the largest sub- 
graph whose members have at least as many neighbors 
as their threshold value ki. This may include finite clus- 
ters as well as any giant component. If the ki are all equal 
we recover the standard fc-core. We define a simple rep- 
resentative example of the heterogeneous-fc-core (HKC) 
in which vertices have a threshold of either 1 or fc > 3, 
distributed randomly through the network with proba- 
bilities / and (1 — /) respectively. This can be directly 
contrasted with bootstrap percolation (BPC), in which 
vertices can be of two types: with probability / they are 
'seed' vertices which are always active, while with proba- 
bility (1 — /) vertices become active only if their number 
of active neighbors reaches a threshold fc. The differ- 
ence between these two processes arises because boot- 
strap percolation is an activation process, beginning from 
a sparsely activated networkjwhile the heterogeneous- 
fc-core is a pruning process [1^ Il5j . beginning from a 
complete graph. It is thus possible to think of these 
two processes as two branches of a hysteresis loop in an 
activation-pruning process. 

We observe two transitions in the size of the giant 
heterogeneous-fc-core (giant- HKC): a continuous transi- 
tion similar to that found in ordinary percolation, and 
a discontinuous, hybrid, transition, similar to that found 
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for the ordinary fc-core. We find a complex phase dia- 
gram for this giant-HKC with respect to the proportion 
of each threshold and the amount of damage to the net- 
work, in which, depending on the parameter region, ei- 
ther transition may occur first. Two similar transitions 
are observed in the phase diagram of the giant component 
of active vertices in bootstrap percolation (giant-BPC). 

Finally, we show that network heterogeneity plays an 
important role. When the second moment of the degree 
distribution is finite but the third moment diverges, the 
giant-HKC (or giant-BPC) appears at a finite threshold 
but not linearly, instead being a higher order transition. 
When the second moment of the degree distribution di- 
verges - as in scale-free networks - the thresholds may 
disappear completely, so that the giant-HKC (or giant- 
BPC) appears discontinuously at a finite value for any 
/ > or p > 0. 



THE HETEROGENEOUS-fc-CORE AND 
BOOTSTRAP PERCOLATION 



Consider an arbitrary, uncorrelated, sparse complex 
network, defined by its degree distribution P{q). In the 
infinite size limit, such networks are locally tree-like, a 
property which enables the analysis we will use. The 
network may be damaged to some extent by the removal 
of vertices uniformly at random. The fraction of surviv- 
ing vertices is p. 

In the heterogeneous-A:-core, each vertex of a network 
is assigned a variable ki £ {0,1,2,...}. The ki values 
are assumed to be uncorrelated, selected from a distribu- 
tion Qk{r)- The heterogeneous-fc-core is then the largest 
subgraph of the network for which each vertex i has at 
least ki neighbors within the heterogeneous-fc-core. To 
find the heterogeneous-Zc-core of a given network, we start 
with the full network, and prune any vertices whose de- 
gree is less than its value of ki. As a result of this prun- 
ing, other vertices will lose neighbors, and may thus drop 
below their threshold, so we repeat the pruning until a 
stationary state is reached. The remaining sub-graph is 
the heterogeneous-fc-core. If it occupies a non-vanishing 
fraction of the original network in the limit that the 
size of the network goes to infinity, we say it is a giant 
heterogeneous-A:-core (giant-HKC) . 

If all ki = 1, then the HKC is simply the con- 
nected component of the network, and the giant-HKC 
is the giant connected component, ex actly as in ordi- 
nary percolation. As is well known l45l |46| this ap- 
pears with a continuous transition at the critical point 
Pc = {q)/[{q^) - {q)], where Y.^<fP{q)■ K k = 2, 

we again have a continuous transition, similar to ordinary 
percolation. If all ki are equal to fc > 3, then we have 
the ordinary /c-core. In this case the giant fc-core appears 
with a discontinuous hybrid transition [^, [TtI |47| . 

Let us briefly discuss the nature of hybrid phase tran- 
sitions. These transitions form a specific new kind of 
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FIG. 1: Top: relative size <Sk of the heterogeneous-fc-core 
(solid curves), which is the subgraph including all vertices 
which meet the threshold requirements of Eq. ([l]), and frac- 
tion <Sb of active vertices in bootstrap percolation (dashed 
curves) as a function of / for the same network - an Erdos- 
Renyi graph of mean degree 5 - with the same fc = 3, at three 
different values of p, corresponding to different regions of the 
phase diagrams Fig. [D 1) p = 0.5, which is between pc and 
Ps for both models. 2) p = 0.61, which is above Ps-k but still 
below ps-b. 3) p = 0.91, which is above pf_k and ps-b. Each 
numbered pair forms the two branches of a hysteresis process, 
followed in the direction marked by the arrows. Bottom: Size 
5gc-k of the giant heterogeneous-fc-core (solid) and size iSgc-b 
of the giant-BPC (dashed) as a function of / for the same 
network and the same values of p. 



phase transition which combines a discontinuity like a 
first order phase transition with a critical singularity 
like a continuous phase transition. In thermodynam- 
ics, where the changes of a control parameter (e.g., de- 
creasing/increasing temperature) are assumed to be in- 
finitely slow, a first order transition has no hysteresis. 
In reality, the changes always occur with a small but 
finite rate, and the 'heating' and 'cooling' branches of 
a first order transition do not coincide, which indicates 
the presence of a metastable state and hysteresis. The 
width of the hysteresis increases with this rate until 
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some limit, which corresponds to the limiting metastable 
state. The resulting limiting curve for the order parame- 
ter has a square- root singularity at the breakdown point. 
For example, if we heat a ferromagnet with a first or- 
der phase transition sufficiently rapidly, the curve 'mag- 
netization M versus temperature T' has a singularity 
M{T) - M{Tb - 0) oc y/Tb - T at the breakdown point 
Tb and discontinuity, that is M {Tb - 0) > M (Tfc 0) = 0. 
After the breakdown, there is no singularity. The suscep- 
tibility also shows a singularity x(T < Tb) oc l/v^Tb^T^ 
and has no singularity above Tb- In this sense, the hybrid 
(mixed) transition is a limiting metastable state for a first 
order phase transition. The important property is that 
the transition is asymmetrical. There are critical fluc- 
tuations and a divergent correlation length on only one 
side of the critical point. Continuous phase transitions 
demonstrate critical fluctuations and divergent correla- 
tion length on both sides of the transition. In first order 
phase transitions there are no critical fiuctuations and 
correlation length is finite everywhere, including at the 
critical point. 

In general we might expect a combination of continu- 
ous and hybrid transitions. To this end we first consider 
the simple case in which the ki are distributed between 
two values, controlled by a parameter /. Specifically: 



Because bootstrap percolation is an activation process, 
while the heterogeneous-/c-core is found by pruning, we 
can characterize them as two branches of the same pro- 
cess, with the difference between the curves shown in 
Fig. [T] indicating hysteresis. Consider beginning from a 
completely inactive network (that may be damaged so 
that some fraction p of vertices remain) . As we gradually 
increase / from zero , under the bootstrap percolation 
process, more and more vertices become active (always 
reaching equilibrium before further increases of /) until 
at a certain threshold value, fci-b a giant active compo- 
nent appears. As we increase / further, the size of the 
giant-BPC traces the dashed curves shown in Fig. [TJ See 
also [4l| . The direction of this process is indicated by the 
arrows on these curves. Finally at / = 1 all undamaged 
vertices are active. Now we reverse the process, begin- 
ning with a fully active network, and gradually reducing 
/, de- activating (equivalent to pruning) vertices that fall 
below their threshold ki under the heterogeneous-fc-core 
process. As / decreases, the solid curves in Fig. [T] will be 
followed, in the direction indicated by the arrows. Notice 
that the size of the giant-HKC for given values of / and p 
is always larger than the giant-BPC. The explanation for 
this difference will be explored in the following Section. 
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for some integer fc > 3. This parameterized HKC has 
as its two limits ordinary percolation (/ = 1) and the 
original fc-core (/ = 0). 

We now contrast this model with bootstrap percola- 
tion. In bootstrap percolation, with probability /, a ver- 
tex is a 'seed' and is initially active, and remains active. 
The remaining vertices (a fraction 1 — /) become active 
if their number of active neighbors reaches or exceeds a 
threshold value fc. Once activated, a vertex remains ac- 
tive. The activation of vertices may mean that new ver- 
tices now meet the threshold criterion, and hence become 
active. This activation process continues iteratively un- 
til a stationary state is reached. The seed and activated 
vertices in bootstrap percolation are analogous to the 
threshold 1 and threshold fc groups in the heterogeneous- 
fc-core. We might expect then that the subgraph formed 
by the active vertices in the stationary state of bootstrap 
percolation might be related to the heterogeneous-fc-core. 
In fact, the two subgraphs are necessarily different, as we 
will describe in detail in Sec. |lTl Nevertheless, the two 
processes have some similar or analogous critical behav- 
iors. Note that we are using a lower threshold of 1 for 
the generalized fc-core, meaning isolated vertices are not 
counted as part of the HKC. We could also use a lower 
threshold of 0, which would include more vertices in the 
HKC, but would yield an identical giant-HKC. For this 
reason, we can compare with bootstrap percolation, in 
which the seed vertices have effectively a threshold of 0. 



Let 5k be the fraction of vertices that are in the 
heterogeneous-fc-core. That is, the total of all compo- 
nents, whether finite or infinite, that meet the threshold 
conditions. This is equal to the probability that an ar- 
bitrarily chosen vertex of the original network is in the 
heterogeneous-fc-core. Let 5gc-k be the relative size of 
the giant heterogeneous-fc-core (that is, the subset of the 
heterogeneous-fc-core which forms a giant component) - 
also the probability that an arbitrarily chosen vertex is in 
the giant heterogeneous-fc-core. The fraction of vertices 
forming finite clusters is therefore iSk — Sgc-k- Note that 
in the standard fc-core this is negligibly small. Similarly, 
let 5b be the fraction of active vertices in the bootstrap 
percolation model, and Sgc-h be the size of the giant com- 
ponent of active vertices. We construct self-consistency 
equations for iSk and iSgc-k in Appendix [X] In networks 
without heavy-tailed degree distributions, that is, whose 
second moments do not diverge in the infinite size limit, 
we find two different transitions for each process. 

We first briefly describe the transitions observed in 
bootstrap percolation, before comparing these with those 
found in the new heterogeneous-fc-core process. For boot- 
strap percolation, above a certain value of p, the giant 
active component (giant-BPC) may appear continuously 
from zero at a finite value of /, /ci-6, and grow smoothly 
with /, sec the dashed lines 1 and 2 in the top panel of 
Fig. [T] For larger p, after the giant active component ap- 
pears, there may also be a second discontinuous hybrid 
phase transition, at fc2-b, as seen in the dashed line 3 
of Fig. [T] There is a jump in the size of the giant ac- 
tive component iSgc-b from the value at the critical point 
(marked by a circle on dashed line 3). When approaching 
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FIG. 2: Left: Phase diagrams for heterogeneous-fc-core (left) in the f-p plane. The giant-HKC is present in regions II and 
III, appearing continuously at the threshold marked by the thin black curve. The hybrid, discontinuous transition occurs at 
the points marked by the heavy black line, beginning from Ps-fc. Above pt-k the first appearance of the giant-HKC is with 
a discontinuous transition. Above Pf-k the giant-HKC appears discontinuously for any / > 0. In region III the vertices with 
threshold k in the HKC form a giant connected component, which appears with a continuous transition (dashed curve). Right: 
Phase diagram for bootstrap percolation. The giant-BPC is present in regions II and HI. The continuous appearance of the 
giant-BPC is marked by the thin solid curve, and the hybrid transition (beginning at Ps-t) by a heavy solid curve. In region 
HI the active vertices with threshold k form a giant connected component. These diagrams are for a Bethe lattice with degree 
5, and for fc = 3, but the diagram for any network with finite second moment of the degree distribution will be qualitatively the 
same. Note that the locations of the continuous transitions from region I to II and from II to HI for bootstrap percolation with 
fc — >■ fc — 1 coincide (up to the point where the discontinuous transition is encountered) with those for the heterogeneous-fc-core. 
The special critical point Ps-t with k reduced by 1 also coincides with ps-k (see also Appendix IX)) . 



from below, the difference of 5gc-b from the critical value 
goes as the square root of the distance from the critical 
point: 

Sbif) = 5b(/c2-b) - a(./c2-b - (2) 

See [4l[ for a complete description of the bootstrap per- 
colation results. 

For the heterogeneous-/c-core, we see a different but 
analogous pair of transitions. Again, for a given p, the 
giant heterogeneous-/c-core (giant-HKC) appears contin- 
uously from zero above some critical value of /, fci-k- 
There may also be a second transition, at fc2-k, where 
again we see a discontinuity in the size of the giant-HKC. 
Now however, the square root scaling occurs as the crit- 
ical point is approached from above: 

5k(/) = 5k(/c2-k) + a{f - /c2-k)'/'. (3) 

See solid line 2 in Fig. [TJ Another important difference 
is that, while for bootstrap percolation the discontinu- 
ous transition is always above the first appearance of the 
giant-BPC - fc2-b > fci-b, but for the heterogeneous- fc- 
core, fc2-k rnay be greater than fci-k - so that the first 
appearance of the giant-HKC is similar to that found in 
ordinary percolation - or less than /ci-fc, with the giant- 
HKC appearing discontinuously from zero, as in the or- 
dinary fc-core BE SI (see Fig. [ID. 



The overall behavior with respect to the parameters 
p and / of each model is summarized by the phase dia- 
grams in Fig. O The diagrams are qualitatively the same 
for any degree distribution with finite second moment. 
Considering first the heterogeneous-fc-core, a giant-HKC 
is absent in the region labelled I. For p below the perco- 
lation threshold Pc, the giant-HKC never appears for any 
/. Above Pc, the giant-HKC appears with a continuous 
transition, growing linearly with / (or p for that matter) 
close to the critical point. The threshold is indicated by 
the thin black line in the Figure, which divides regions 
I and II. Compare line 1 of Fig. [TJ From Ps^k a second, 
discontinuous, hybrid transition appears. This is marked 
by the heavy black curve in Fig. [21 As already noted in 
Eq. ^ the size of the giant-HKC grows as the square 
root of the distance above this second critical point. At 
the special point Ps-k, however, the size of the discon- 
tinuity reduces to zero, and the scaling near the critical 
point is cube root: 

5k(/) = 5k(/c2-k) + a{f - /c2-k)'/^ (4) 

At first the hybrid transition occurs after the continuous 
appearance of the giant-HKC (see line 2 of Fig. [1]), but at 
Pt-k the two transitions cross, and the giant-HKC begins 
to appear immediately with a jump. Finally at Pf^k, the 
giant-HKC begins to appear discontinuously immediately 
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from / = (line 3 of Fig. [T]). Thus, the first appearance 
of the giant-HKC has a classical percolation like tran- 
sition below pt-k, while above pt-k the appearance is 
similar to that found in the ordinary k-coie. Within the 
heterogeneous- /c-core, the vertices with threshold k may 
form a giant component themselves. This occurs in the 
region labelled III. This giant component appears with a 
continuous transition. 

For bootstrap percolation, the giant active component 
again appears only above the percolation threshold Pc- It 
also appears at first with a linear, continuous transition 
(when the degree distribution has finite third moment), 
but at a larger value of /. The point of appearance is 
marked by the thin curve in Fig. [2] which divides regions 
I and II. Again, a second transition appears for larger p, 
beginning at the special critical point Ps-t- This tran- 
sition is marked by the heavy solid curve in the Figure. 
See also line 3 of Fig. [TJ The scaling near the hybrid 
transition is again square root, but this time only when 
approaching the transition from below. At the special 
critical point Ps-b, the scaling becomes cube root: 

SM) = '5b(/c2-b) - a(/c2-b - (5) 

Note the difference between Eqs. ^ and ([3]) and between 
([5]) and In bootstrap percolation, the hybrid transi- 
tion always occurs above the continuous one, and neither 
reaches / = 0. Again, the active threshold k vertices 
form a giant component in Region III. Note also that the 
special critical point Ps-t > Pf~k, so that for a givenp we 
may have a hybrid transition for the HKC or for BPC, 
but not for both. It turns out that the special critical 
point Ps-b for bootstrap percolation whose value of k is 
one less coincides with the special critical point for the 
heterogeneous- fc-core. Furthermore, the location of the 
continuous appearance of the giant heterogeneous-fc-core 
for a given k also coincides with the appearance of the 
giant component of bootstrap percolation for fc — 1. The 
same is also true for the appearance of the giant com- 
ponents of vertices with threshold k. This is clear from 
the equations given in Appendix [X] If the value of k is 
increased, the locations of the hybrid transitions move 
toward larger values of p, and there is a limiting value 
of k after which these transitions disappear altogether. 
For both the heterogeneous-fc-core and bootstrap perco- 
lation, this limit is proportional to the mean degree. Note 
also that the continuous transition also moves slightly 
with increasing fc, and in the limit fc — cx), tends to the 
line pf = Pc for both processes. 

II. SUBCRITICAL CLUSTERS, CORONA 
CLUSTERS, AVALANCHES 

It is clear from Figs. [T] and [5] that even though boot- 
strap percolation and the heterogeneous-fc-core described 
above have the same thresholds and proportions of each 
kind of vertex, the equilibrium size of the respective gi- 
ant components is very different. The difference results 



from the top-down vs bottom-up ways in which they are 
constructed. To find the heterogeneous-fc-core, we be- 
gin with the full network, and prune vertices which don't 
meet the criteria, until we reach equilibrium. In contrast, 
bootstrap percolation begins with a largely inactive net- 
work, and successively activates vertices until equilibrium 
is reached. To see the effect of this difference, we now de- 
scribe an important concept: the subcritical clusters of 
bootstrap percolation. 




FIG. 3: Top row: Subcritical clusters of different sizes in boot- 
strap percolation. Left: Because they start in an inactive 
state, two connected vertices (shaded area) cannot become 
active if each has k — 1 active neighbors. (In this example 
A: = 3.) The same follows for clusters of three (center) or 
more vertices (right). If any member of a subcritical cluster 
gains another active neighbor, an avalanche of activations en- 
compasses the whole cluster. Bottom row: Similar clusters 
would be included in the heterogeneous-fc-core. 




FIG. 4: Corona clusters of different sizes in heterogeneous- 
fc-core. Left: Because they are included unless pruned, two 
connected vertices (whose threshold is fc) form part of the 
heterogeneous-fc-core if each has fc — 1 other neighbors in the 
core, as each is 'assisted' by the other. In general (right) 
a corona cluster consists of vertices (with threshold fc) that 
each have exactly fc active neighbors, either inside or outside 
the cluster. If one neighbor of any of the cluster vertices is 
removed from the core (for example, the one indicated by the 
arrow), an avalanche is caused as the entire cluster is pruned. 



A subcritical cluster in bootstrap percolation is a clus- 
ter of activatable vertices (i.e. not seed vertices) who 
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each have exactly k — 1 active neighbors external to the 
cluster. Under the rules of bootstrap percolation, such 
clusters cannot become activated. The vertices within 
the cluster block each other from becoming active - see 
Fig. [3l Now compare the situation for the heterogeneous- 
k-coie. Any cluster of threshold k vertices which each 
have k — 1 neighbors in the core external to the clus- 
ter, is always included in the heterogeneous-fc-core. for 
the heterogeneous-fc-core, vertices in clusters like those 
in Fig. m assist one another. Thus the exclusion of sub- 
critical clusters from activation in bootstrap percolation 
accounts for the difference in sizes of the bootstrap per- 
colation core and the heterogeneous-fc-core. 

Subcritical clusters have another important property, 
that helps us to understand the discontinuity at the sec- 
ond transition. A vertex in a subcritical cluster becomes 
active if it gains an extra active neighbor (for example 
through an infinitesimal change in p or /). This in turn 
allows each of its neighbors in the cluster to activate. A 
domino-like effect ensues, leading to an avalanche of ac- 
tivations - one extra active neighbor of any vertex in the 
subcritical cluster leads to the whole cluster becoming 
active. Thus the rate of change of Sh is related to the 
sizes of the subcritical clusters. Almost everywhere, if we 
choose a subcritical vertex at random, the mean size of 
the subcritical cluster to which it belongs is finite. How- 
ever, exactly at the second threshold, the mean size of 
the subcritical cluster to which a randomly chosen ver- 
tex belongs diverges as we approach from below. This 
was shown in [4lj. Thus, approaching this point, an in- 
finitesimal decrease in / (or p) leads to a finite fraction of 
the network becoming activated, hence a discontinuity in 
5b (and also in 5gc-b)- The distribution of avalanche sizes 
near the transition is determined by the size distribution 
G{s) of subcritical clusters, which at the critical point 
can be shown to follow G{s) ~ s~^^^. This can be shown 
using a generating function approach, as demonstrated 
in [m . A similar method can be found in [13, 147l - l49| . 

We can then understand the hybrid transition in the 
heterogeneous-fc-core, by considering the relevant clus- 
ters with similar properties, the corona clusters. Corona 
clusters are clusters of vertices with threshold fc that 
have exactly fc neighbors in the HKC. These clusters are 
part of the heterogeneous-fc-core, but if any member of a 
cluster loses a neighbor, a domino-like effect leads to an 
avalanche as the entire cluster is removed from the HKC. 
The corona clusters are finite everywhere except at the 
discontinuous transition, where the mean size of corona 
cluster to which a randomly chosen vertex belongs di- 
verges as we approach from above [1, [13, H^]. Thus, an 
infinitesimal change in / (or p) leads to a finite fraction 
of the network being removed from the heterogeneous-fc- 
core, hence a discontinuity in 5k (and also in 5gc-k)- The 
size distribution of corona clusters, and hence avalanches 
at this transition also goes as a power law with exponent 
-3/2. 



III. SCALE-FREE NETWORKS 

The results above, in Sec. U and in Figs. [1] and [2l are 
qualitatively the same for networks with any degree dis- 
tribution which has finite second and third moments. 
When only the second moment is finite, the phase dia- 
gram remains qualitatively the same, but the critical be- 
havior is changed. Instead of a second order continuous 
transition, we have a transition of higher order. When 
the second moment diverges, we have quite different be- 
havior. 

To examine the behavior when the second and third 
moments diverge, we consider scale-free networks, with 
degree distributions tending to the form 

P(g) « q-'' (6) 

for large q. At present we consider only values of 7 > 2. 

To find the behavior near the critical points, we begin 
with a self-consistency equation for X , which is the prob- 
ability that an arbitrarily chosen edge leads to an infinite 
{ki — l)-ary tree (see Appendix El) The equation is given 
in the Appendix, as Eq. (IA4[) . The probability iSgc-k can 
be written in terms of X, and in fact both X and iSgc-k 
grow with the same exponent near the appearance of the 
giant active component. We expand the right hand side 
of Eq. (IA4[) near the appearance of the giant-HKC (that 
is, near X = 0). When 7 < 4, the third and possibly 
second moment of the degree distribution diverge. This 
means that coefficients of integral powers of X diverge, 
and we must instead find leading non-integral powers of 
X. 

When 7 > 4, the second and third moments of the 
distribution are finite, so the behavior is the same as 
already described in Sec. HI Eq. (jA4p leads to 

X = ciX + C2X^ + higher order terms , (7) 

which gives the critical behavior X gc {f — fc2-k)^ with 
13 = 1. When 3 < 7 < 4, the linear term in the expansion 
of Eq. (|A4[) survives, but the second leading power is 
7-2: 

X = ciX + C2X'^-^ + higher terms , (8) 

where the coefficients ci and C2 depend on the degree 
distribution, the parameters p and /, and the (non-zero) 
value of Z. The presence of the linear term means the 
giant-HKC appears at a finite threshold, but because the 
second leading power is not 2, the giant-HKC grows not 
linearly but with exponent /? = 1/(7 — 3). This means 
that the phase diagram remains qualitatively the same 
as Fig. [21 however, the size of the giant-HKC grows as 
(/ ~ /ci)^ with /3 = 1/(7 — 3). This is the same scaling 
as was found for ordinary percolation [5C| . 

For values of 7 below 3, the change in behavior is more 
dramatic. When 2 < 7 < 3, the second moment of P{q) 
also diverges, meaning the leading order in the equation 
for X is no longer linear but 7 — 2: 

X = diX'>^^ + higher terms . (9) 
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From this equation it follows that there is no threshold 
for the appearance of the giant- HKC (or giant-BPC). The 
giant-HKC appears immediately and discontinuously for 
any / > (or p > 0), and there is also no upper limit 
to the threshold k. This behavior is the same for boot- 
strap percolation, so the (featureless) phase diagram is 
the same for both processes, even though the sizes of the 
giant-HKC and giant-BPC are different. 



IV. DISCUSSION 

We have introduced a new concept, the heterogeneous- 
k-coie, an extension of the well known fc-core of com- 
plex networks. A simple representative example of the 
heterogeneous- fc-core (HKC) has vertices with randomly 
assigned thresholds of either 1 (with probability /) or 
fc > 3 (with probability 1 — /). This heterogeneous- fc- 
core has a complex phase diagram, including two types 
of phase transition: a continuous transition at the ap- 
pearance of the giant heterogeneous- fc-core, similar to the 
ordinary percolation transition, and a second, discontin- 
uous, hybrid phase transition. This second transition 
is similar to that found for the ordinary fc-core, but it 
may occur after the first continuous appearance of the 
giant-HKC or before. The first transition occurs when 
the vertices of both kinds reaching their threshold form 
a giant percolating cluster. The second transition occurs 
when the mean size of avalanches of pruned vertices di- 
verges. This can be understood by considering corona 
clusters (clusters of vertices which exactly meet the up- 
per threshold). The mean size of the corona cluster to 
which an arbitrarily chosen vertex belongs diverges as we 
approach the second transition from above. The size of 
pruning avalanches are determined by these corona clus- 
ters, and so the transition is discontinuous. 

We have contrasted the heterogeneous-fc-core with 
bootstrap percolation, in which there are also two kinds 
of vertices, but the core is defined by an activation pro- 
cess, rather than a pruning process. The phase dia- 
gram for bootstrap percolation therefore does not coin- 
cide with that of the heterogeneous-fc-core. Furthermore, 
the two processes can be thought of as two branches 
of an activation-pruning process, forming a hysteresis 
loop. The difference between the giant heterogeneous- 
fc-core and the giant bootstrap percolation component 
results from the subcritical clusters of bootstrap percola- 
tion. These clusters cannot be activated as all members 
have fc — 1 active neighbors outside the cluster. However, 
the equivalent clusters would be included in the giant 
heterogeneous-fc-core. 

All of these results are strongly dependent on network 
structure. If the third moment of the degree distribution 
is finite, we obtain the results just described, with the 
giant-HKC growing linearly above the continuous thresh- 
old. This is the case if the degree distribution decays 
faster than a power-law q~'^ with 7 > 4 for large degree 
q. If instead the degree distribution tends to a power-law 



with 3 < 7 < 4, so that the second moment is finite while 
the third moment diverges, the phase diagram is qualita- 
tively the same, but the continuous transition is of higher 
order. If 2 < 7 < 3, the second moment of the degree 
distribution diverges, and the situation is more extreme. 
The giant-HKC appears immediately at a finite size for 
any / > or p > 0, showing that, in common with be- 
havior found in other systems, such scale-free networks 
are extremely resilient to damage. 



Acknowledgments 

This work was partially supported by the follow- 
ing projects PTDC: FIS/71551/2006, FIS/108476/2008, 
SAU-NEU/103904/2008, and MAT/114515/2009, and 
also by the SOCIALNETS EU project. 



Appendix A: Self-Consistency Equations 

Here we construct the self-consistency equations that 
the probabilities 5k, 5gc-k, ^b, and 5gc-b must obey. The 
(usually numerical) solution of these equations lead to the 
phase diagrams and other results presented in Section ID 
We have already given them for the case of bootstrap 
percolation [4l| . 

• o 

threshold 1 threshold k 




O © ^gc-k 

TABLE I: Symbols used in graphical representations of self- 
consistency equations for the heterogeneous-fc-core. 

The probability iSk that an arbitrarily chosen vertex 
belongs to the heterogeneous-fc-core (HKC) is the sum 
of the probabilities that it has ki — 1 and at least one 
neighbor in the core, or fci = fc and has at least k neigh- 
bors in the core. We can represent this diagrammatically 
as: 

We define Z in terms of a '(fc^ — l)-ary tree', a gener- 
alization of the (fc — l)-ary tree. A (fc^ — l)-ary tree is 
a sub-tree in which, as we traverse the tree, each vertex 



8 



encountered has at least ki — 1 child edges (edges leading 
from the vertex, not including the one we entered by). 
The variable ki can be different at each vertex. In our 
example, ki is either 1, in which case the vertex does not 
need to have any children (though it may have them) , and 
is only required to be connected to the tree, or ki = fc, in 
which case it must have at least fc — 1 children. 

The probability Z can then be very simply stated as 
the probability that, on following an arbitrarily chosen 
edge in the network, we reach a vertex that is a root of 
a {ki — l)-ary tree. For the specific case considered in 
this paper, the vertex encountered either has fci = 1, or 
it has fc — 1 children leading to the roots of {ki — l)-ary 
trees. The probability Z is represented by a square in 
the diagram. A bar represents the probability {1 — Z), a 
black circle represents a vertex with fc^ = 1, and a white 
circle a vertex with ki = k - see TableHl These conditions 
can be written as binomial terms, and summing over all 
possible values of the degree of i, this diagram can be 
written in mathematical form as: 



s^=pfY,p{q)Y,r\z\i-zy 

g=l 1=1 ^ ^ 



similar to Z , but with the extra condition that the sub- 
tree reached must extend indefinitely. We represented X 
by an infinity symbol (Table The diagram for iSgc-k 
is: 




which is equivalent to the equation: 

oo / \ 

q=0 m=l ^ ^ 



+p{i-f)j2P{q)j2(,){i-zr-' 

q—k l—k 



tQx"(z-xf 

m=l ^ ' 



(A3) 



To find X, we construct a self-consistency equation 
from the diagram 



p(i - /) ^ p{q) E ( ; ) ^' (1 - ^y^' ' (Ai) 



where the factor p accounts for the probability that the 
vertex has not been damaged. 

To calculate Z, we construct a recursive (self consis- 
tency) expression in a similar way, based on the definition 
given above. This is represented by the diagram: 




leading to 



>1 \J^-'~ 



n=0 m=l ^ ^ 




>>k- 1 



which, in equation form is: 



^ (q) ^ V ^ 

q>k l = k-l ^ 



-^^{Z,pJ) 



(A2) 



We have used that qP{q) / (q) is the probability that the 
vertex reached along an arbitrary edge has degree q. 
Solving Eq. (|A2|) (usually numerically) for Z and then 
substituting into Eq. ()Aip allows the calculation of 5k. 

We follow a similar procedure to calculate the size 
of the giant heterogeneous-fc-core (giant- HKC), which is 
equal to the probability Sgc-k that an arbitrarily chosen 
vertex is a member of a heterogeneous-fc-core component 
of infinite size. We denote by X the probability that an 
arbitrarily chosen edge leads to a vertex which is the root 
of an infinite (fc.; — l)-ary tree. That is, the definition is 



m— 1 



(A4) 



Solution of Eqs. (jA2l) and (IA4[) then allows the calcula- 
tion of Sgc-k through Eq. (jA3p . Note that when there are 
multiple solutions of Z, we choose the largest solution as 
the 'physical' one. 

To find the appearance of the giant component for a 
given p, we find leading terms for X 1, and solve for 
/, as described in Sec. IIIII To calculate the location of 
the hybrid transition we note that at this critical point a 
second solution to Eq. (jA2p appears. This occurs when 
the function "^{Z) just touches the line Z, which must 
be at a local extremum of ^f/Z: 



(A5) 



Expanding Eq. (jA2[) about Zc, the value of Z at the crit- 
ical point (at the top of the jump), and using Eq. (jA5l) . 
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we see that Z grows as the square-root of the distance 
from the critical point. Using Eq. (|A1|) we find Eq. ([3]). 

Furthermore, at the special point Ps-k where the sec- 
ond transition disappears, by a similar argument, a fur- 
ther condition must also be satisfied: 



dZ2 I Z 



0. 



(A6) 



Thus, the critical point p^-fc is determined by simulta- 
neous solution of Eqs. (|X2|) . (IX5|) . and (IM)) . This in 
turn leads to cube root scaling above the threshold, hence 
Eq. Q. 

For bootstrap percolation, we can construct similar 
self-consistency equations in order to calculate 5b and 
^gc-b- Let Y be the probability (counterpart of Z) that 
on following an arbitrary edge, we encounter a vertex 
that is either a seed or has k active children. As dis- 
cussed in Section |TT1 activation in bootstrap percolation 
must spread through the network, meaning that the ver- 
tex needs fc active downstream neighbors in order to be- 
come active (and thus provide an active neighbor to its 
upstream 'parent'). Repeating the diagrammatic method 
described above, we arrive at the equation: 



(q) ^ V I 

q>k+l l = k 



(A7) 



Note that Eq. (IA7[) differs from (jA2l) because the number 
of children required is k not fc — 1. This is equivalent 
to excluding the subcritical clusters represented in Fig. 
131 As an aside, consider the probability T^sub that an 
arbitrary edge leads to a vertex in a subcritical cluster. 
This is clearly simply the probability that the vertex has 
exactly k — 1 active neighbors, thus 

= p(i - /) E ^ : ;) r-Hi - 

(A8) 

Comparing Eqs. (|A2[) and (|A7p . we see that the right 
hand side of Eq. (jA8l) contains precisely the terms that 
are counted in Eq. (|A2|) but absent from (|A7p . Thus 
clusters of vertices all having exactly k — 1 active neigh- 
bors are excluded from the active component in boot- 
strap percolation, while vertices having k — 1 neighbors 
in the heterogeneous-fc-core are always included in the 
heterogeneous-fc-core. Of course because of the self- 
recursion, the value of Z is necessarily different from that 
of Y by more than just the amount of these terms. 

Drawing diagrams similar to those given for the 
heterogeneous-Zc-core allow the construction of further 
self-consistency equations for the remaining quantities of 
interest. The probability that an arbitrarily chosen ver- 
tex is active, Sh is then identical to Eq. (jAip but with Y 



replacing Z: 

oo q 

5b-p/E^(9)E(;)^'(l-^)' 



9=1 



1 = 1 



q=k 



l=k 



p(l-/)^P(g)^(^)r'(l-y)''-'. (A9) 



The equation for iSgc-b follows similarly. We introduce the 
probability W that, upon following an arbitrary edge, we 
reach a vertex that is active and also has at least one edge 
leading to an infinite active subtree. Then W obeys: 

q=0 m=l ^ ^ 



l=k 



I 



J2 (^^w"'{Y-wy-"' 



(AlO) 



which again, differs from Eq. (|A4[) in that the limit is 
k — 1 not k. Then 5gc-b obeys: 

oo ' / \ 

g=0 m=l ^ ^ 



q—k l—k 



J2 jW^'iY-wy-' 



(All) 



which is identical in form to Eq. (|A3I) , though the values 
of Y and W (for bootstrap percolation) will be different 
from those of Z and X (for the heterogeneous- fc-core). 
Note also that, for bootstrap percolation the physical so- 
lution for Y is always the smallest of Eq. (|A7p . The 
discontinuous transition occurs at the point where: 



5F y'=° 



(A12) 



and at the special point Ps-b there is one more condition. 



dY^ V Y 



T7 =0, 



(A13) 



where the function ^(Y) is defined by Eq. (jA7[) . 

Finally we note that the appearance of the giant com- 
ponent of threshold k vertices in the heterogeneous- fc- 
core process can be found in a similar way to the appear- 
ance of the giant- HKC. We define R to be the probability 
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that an arbitrarily chosen edge leads to the root of an in- 
finite subtree that is a (fc^ — l)-ary tree with all of the 
fcj = 1 vertices removed. This probability then obeys a 
self consistency equation similar to that for X : 

q=k /=fe-l ^ ^ 



t(„>"(z-«)' 

m=l ^ ^ 



(A14) 



The appearance of the giant component of threshold k 
vertices is then found by expanding this equation to lead- 
ing order with respect to R and solving for /. This leads 
to an equation which gives the dashed line in Fig. [21 Re- 
call that Z{f,p) is determined by Eq. (jA2p . A similar 
procedure yields the corresponding transition in boot- 
strap percolation. 



Appendix B: General Form of Equations 

In this paper, we have examined only a special case of 
the heterogeneous-fc-core, in which vertices have thresh- 
old cither 1 or fc > 3. For completeness, we now give the 
self-consistency equations for arbitrary threshold distri- 
bution Q{r). The size Sk of the heterogeneous-fc-core is 



5k=p5]Q(r)^P(g) 



q—r 



.l^r 



q-l 



(Bl) 



where, as above, Z is the probability of encountering a 
vertex i which is the root of a (fc^ — l)-ary tree: 



r>l q=r l=r-l 



9-1 



9-1 



z\i-zy-^-K 

(B2) 



Similarly, the equation for the size of the giant active 
component is 



r>l q=r l=r ^ ^ 

X E {DlX"\Z-Xr^\ (B3) 



m— 1 



where X obeys: 



9-1 



(9) 



l=r-l 



(1 - zy 



-i-i 



I 



(B4) 



We do not derive any results for this general case, but 
we can speculate that a more complicated phase diagram 
would appear. If any vertices have threshold less than 
3, i.e. Q(l) -|- (5(2) > 0, we would find a continuous 
appearance of the giant-HKC. Thresholds of 3 or more, 
on the other hand, contribute discontinuous transitions, 
and it may be that there are multiple such transitions. 
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